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We examined the bonding and magnetic properties, as well as the stability, of three magnetic 
half-Heusler alloys, namely LiMnZ, with Z=N, P or Si, in the three different atomic ordering of 
the Cli, crystal structure (i.e. a, /3, and 7 phases). At their respective optimized lattice constant, 
each material in the primitive cell is a ferromagnetic metal. Assuming a proper matching substrate 
could be found, a-LiMnN, LiMnSi in the ft and 7 phases, and LiMnP in the a and f3 phases can 
be ferromagnetic half-metals at lattice constants larger than their optimized values. We also found 
that a-Lio.75MnP, /3-Lio.7sMnZ, and pnictide 7-Lio.75MnZ can be half-metals at lattice constants 
slightly larger than the optimized values. Subsequently, we found that ferromagnetic LiMnZ obeys 
a modified Slater-Pauling curve, and all ferromagnetic Lio.7sMnZ, except for in the a phase, obey a 
similar curve. By comparing /3-LiMnP to the metastable zinc blende phase of MnP we determined 
the role of Li in the structure with respect to the elastic stability, electronic properties, and magnetic 
properties. Finally, we examined the possibility of antiferromagnetic ordering in these materials and 
found that by using the tetragonal unit cell, consisting of two formula units, LiMnN and LiMnP 
favor antiferromagnetic alignment. 

PACS numbers: 85.75.-d,75.50.Ee,71.15.Mb 
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I. INTRODUCTION 

Since de Groot et alP predicted that NiMnSb would 
exhibit half-metallic properties, this kind of ternary com- 
pound has attracted much attention owing to the spe- 
cial properties of the half-metallic phase. The formula, 
XYZ, characterizes these ternary compounds which are 
denoted as half-Heusler alloys or semi-Heusler alloys. 2 
The crystal structure of these compounds, Clb, is sim- 
ilar to a full-Heusler alloy (X 2 YZ), or L2i, but missing 
one X. For half-Heusler alloys, there are three different 
atomic arrangements, shown in Fig. [I] used by different 
groups, called a, /3 and 7 phases P In Table [TJ the posi- 
tions occupied by the three atoms and the vacancy are 
labeled according to the notation defined by Wyckoff. 4 
The arrangement determined by Mastronardi et al. P us- 
ing x-ray diffraction, is the 7 phase. In this structure, the 
non-magnetic element Z is nearest neighbor (nn) to the Y 
atom. Mancoff et alP grew superlattices of half-Heusler 
alloys in the a phase. In this case, the Z atom has Y as 
the second neighbor (sn) and the X atom is nn to both 
Z and Y. The a phase is mainly used by de Groot et alP 
and Galanakis et alP Finally, the f3 phase, described by 
Larsen et al.P has Z nn to both X and Y. 

One of the authors (C. F.) examined lithiated half 
Heusler alloys in the (3 phase P Examples are LiMgZ 
(Z=N P, As, Bi), LiZnP, LiCdP and LiAlSi. The co- 
valent bonding between the Y and Z atoms forms the 
gap of the compound. Since there are no magnetic ele- 
ments in any of these alloys, several interesting issues rel- 
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TABLE I. The positions of the three atoms, X, Y and Z, and 
the vacancy in terms of the Wyckoff notation: 4a = (0, 0, 0)a, 
46 = (V2,0,0)a, 4c = (1/4, V«, V 4 ) a an d 4d = (1/4, V 4 . 3 / 4 K 
where a is the lattice parameter. The zinc blende structure, 
YZ, is included for comparison. 



evant to the magnetic properties were not addressed. We 
therefore decided to examine three simple lithiated half- 
Heusler alloys involving the magnetic element, Y=Mn, 
namely LiMnN, LiMnP and LiMnSi. Jungwirth et al. 
incorporated a magnetic element into a Li based half- 
Heusler alloy by epitaxially growing /3-LiMnAs on an 
InAs substrateP Using two formula units, LiMnAs is an- 
tiferromagnetic in the ft phase and thus has no magnetic 
moment. Two formula units of XMnZ are necessary in 
the antiferromagnetic case to allow the Mn to align with 
antiparallel spins. Half-metallic fully compensated (FC) 
fcrrimagnets also have no magnetic moment but are semi- 
conducting as well. Half-metals, on account of their large 
magnetic moments, and FC half-metals, because of their 
lack of stray magnetic fields, are generally anticipated to 
play important roles in spintronic device applications. 




FIG. 1. (Color) Three phases of the half-Heusler alloys, LiMnZ. The sizes of the spheres relate to the covalent bonding radius 
of Li, Mn and Z=P, respectively. 



Using the a structure of transition metal (TM) based 
half-Heusler alloys (X=Fe, Co, Ni; Y=TL V, Zr, Nb, 
Mn), Tobola et alP^ and Galanakis et alP showed that 
the d-d hybridization between X and Y, of half metallic 
CrMnSb, is responsible for the states at the gap in the in- 
sulating channel. They did not concern themselves with 
the d-p bonding states between X and Z. We suggested 
an alternate point of viewP^ Z is the most electronega- 
tive among the three atoms so it should determine the 
bonding properties. From this point of view, the TM 
based half-Heusler alloys, such as CrMnSb, in the a and 
7 structures differ by which of the two TM atoms is the 
nn of Z. Their properties differ accordingly. 

How does our suggestion, to consider the nn of Z, work 
for Li-related half-Heusler alloys with the magnetic ele- 
ment in LiMnZ? In the a structure, Z is nn to Li; in 7, 
is it nn to Mn; and Z is nn to both Li and Mn in the 
P phase. Furthermore, Li does not have any d-states, 
and it easily gives up its outermost electron, so the three 
structures are expected to have a more profound differ- 
ence than the TM half-Heusler alloys. We checked that 
our calculations agreed with this prediction by contrast- 
ing the electronic and magnetic properties of LiMnZ, zinc 
blende (ZB) MnZ, and Li .7sMnZ elements. 

Since the local moment determines the saturation 
magnetization — which in turn relates to the Curie tem- 
perature Tc — large magnetization can lead to high Tc- 
As shown in CrAs, the measured saturation magnetiza- 
tion equates to 3 /^b /unit cell and the measured Tq is 
480 K. 12 It is, therefore, advantageous for a compound 
to have a large magnetic moment. Based on the sim- 
ple ionic modelp^ Mn should transfer three of its valence 
electron to the pnictogens or four electrons to Si. The 
remaining electrons at the Mn site should form a local 
magnetic moment of 4 /xb /unit cell (for the pnictides, N 
and P) or 3^B/unitcell (for LiMnSi). If the Li gives its 
electron to Mn, there is a possibility of increasing the 
local moments at Mn by one. Can these half-Heusler 
alloys have magnetic moments greater than or equal to 
4 /i B /unit cell? 

We also want to determine if these alloys exhibit half- 
metallic properties. A half-metal has one spin channel 



with metallic properties, while the oppositely oriented 
spin channel displays the characteristics of an insulator. 
The criteria for the half metallic phase from a theoreti- 
cal point of view are: (i) the density of states (DOS) at 
the Fermi energy (Ep) should be non-zero in the metallic 
channel and zero in the insulating channel, and (ii) the 
calculated magnetic moment per unit cell (or formula 
unit, f. u.) should be an integer!^ Perfect spin polariza- 
tion at Ep is ideal for devices employing magnetioresis- 
tance. 

The generalized Slater-Pauling curve for TM-based 
Heusler alloys, proposed by Kiiblerf^ has had amazing 
success predicting magnetic moments, especially in fully 
compensated (FC) half-Heusler alloysp^ where there is 
no net magnetic moment. In the FC configuration, differ- 
ent elements align antiparallel with similar magnitudes.^ 
This is distinct from the antiferromagnetic configuration 
where the same elements align antiparallel. The gen- 
eralized Slater-Pauling curve, modified from the Slater- 
Pauling electron counting rulej 18 * 19 ^ governs the magnetic 
moments of 3-d TM metals and their binary alloys by 
considering the results from band structure calculations. 
For TM-based half-Heusler alloys, the curve predicting 
the contribution due to the local moment is 



M=(iV t -18)/x B /f.u. 



(1) 



where N t is the total number of valence electrons per f. u. 
and the number 18 is deduced by the average number of 
down-spin valence electrons per atom = 3 from band 
structure calculations. These electrons pair with up-spin 
electrons and do not contribute to the magnetic moment. 
There are three atoms per f. u., thus the total number of 
paired electrons is 3 x (2 x nj.) = 18. If a TM half- 
Heusler alloy has a total of 18 valence electrons per f. u., 
Eq. predicts no moment and the alloy is FC. The 
first-principle calculations on CrMnSb by §a§iogltP^ and 
Shaughnessy et alP^ agree with this prediction, however, 
we recently reported that there are exceptions when the 
moments predicted by the curve are negative.^ For Li- 
based half-Heusler alloys, Eq. ([I]) must be altered since 
the total number of electrons is different, thereby caus- 
ing a different band structure. A new expression, based 
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on the Li-based half-Heusler alloy band structure calcu- 
lation, will be given and serves as a precursor for future 
predictions. Furthermore, a reliable prediction can facil- 
itate the selection of elements to form new half-Heusler 
alloys, in particular ones with large magnetic moments. 

Finally, for these alloys to be useful in devices, the 
issue of stability should be addressed. Intuitively, the /? 
structure is similar to the ZB structure, except for the 
presence of the Li atom. We compared the stability of 
/3-LiMnP to meta-stable ZB MnP by comparing phonon 
spectra determined by a response function. 

In this paper, we address the following issues: 

A. How do the bonding properties of each arrangement 
differ? 

B. Can the three half-Heusler alloys have large mo- 
ments in either ferromagnetic or half metallic 
phases? 

C. How should the Slater-Pauling curve be altered for 
the lithiated alloys? 

D. Can the three alloys be fully compensated (FC) 
half-metals or antiferromagnets? 

E. How does Li play its role in the structures? 

F. Can the lithiated /3-phase half-Heusler alloys be 
more stable than their ZB counterparts? 

In Section [H] we discuss the models used to answer these 
questions, and present a brief description of the methods 
of calculation. Results and discussion of the above issues 
will be presented in Section |III| Finally, in Section |IV[ 
we summarize our findings. 

II. STRUCTURAL MODELS AND METHOD OF 
CALCULATION 

To address the issues mentioned above, we have to con- 
sider primitive, tetragonal and conventional cubic cells 
for each alloy. The primitive cell consists of the three 
atoms, Li, Mn and Z, with lattice vectors connecting 
the corner (0, 0, 0) to the face centers of the conventional 
cube. We use the primitive cell to find the optimized lat- 
tice constants by means of minimizing the total energy of 
each alloy with respect to the lattice constant. The prim- 
itive cell is also used for the response function phonon 
calculation of /3-LiMnP and MnP. To search for half- 
metals, and to determine the role of Li in the electronic 
and magnetic properties, we used the conventional cubic 
cell. This cell consists of four formula-units for LiMnZ. 
Lio.ysMnZ is identical except missing one Li. In other 
words, we model LiMnZ as Li4Mn4Z4 and Lip^MnZ as 
LisMn4Z4, respectively. Since these alloys have one mag- 
netic transition metal element in the primitive cell, we 
do not expect them to show FC properties. However, 
it is possible to calculate the antiferromagnetic proper- 
ties. We constructed a tetragonal unit cell, consisting of 



two formula- units of LiMnZ, to allow the possibility of 
antiparallel alignment of the same species of atoms. 

We used the spin polarized version of the vasp 
code^HHH] w hich is based on density functional theory^! 
(DFT). The generalized gradient approximation (GGA) 
of Perdew et alP (PBE) was used to treat the exchange- 
correlation between electrons. Except for the value of the 
semiconducting gap, GGA provides realistic bonding and 
magnetic properties. The value of the semiconducting 
gap and the robustness of the half-metallicity are not cru- 
cial issues at this point, so we did not consider first princi- 
pal local density approximation plus dynamic mean field 
theory (LDA+DMFT) or GW methods. These meth- 
ods improve upon the conduction states in semiconduc- 
tors and half-metals,^ and the LDA+DMFT method has 
been used to show that half-metallic properties may dis- 
appear at finite temperaturesPZl 

T he V ASP package provides projector-augmented- 
wave j28i29| ( PAW ) potentials for Li, Mn, N, P and Si that 
were constructed using PBE. We used a basis of plane 
waves with a 1000 eV kinetic en ergy cutoff for all calcu- 
lations. The Monkhorst and PadP^ meshes of (17, 17, 17) 
and (11, 11, 11) were adopted for the Brillouin Zone (BZ) 
of the primitive and the conventional cubic cells, respec- 
tively. The tetragonal unit cell used a mesh of (17, 17, 17) 
as well. Using these values, the convergence of the total 
energy and the magnetic moment of any sample are bet- 
ter than 1.0 meV and l.Om/iB, respectively. 

To address the stability of the lithi ated compounds, 
we used the ABINIT software package^HS] to preform 
response- function phonon calculations. We used the same 
exchange-correlation functional as in vasp and used com- 
parable convergence parameters to calculate the ground 
state structure. To calculate the force constants, a 
4x4x4 irreducible q-point grid, centered at q = (0, 0, 0) 
was constructed. 



III. RESULTS AND DISCUSSION 

A. How do the bonding properties of each 
structure differ? 

Using the primitive unit cell, we first determined the 
optimized lattice constants of the three compounds in 
the three phases. The results are summarized in Ta- 
ble |TTJ The lattice constants correlate with the cova- 
lent radii of Z (r cov (N) = 0.71 A, r cov (P) = 1.09 A, 
and r co „(Si) = 1.14 A) and the total energies are de- 
termined by the electronegativity of Z (x(P) = 2.19, 
X(N) = 3.04, x(Si) = 1.90 compared to x(Mn) = 1.55 
and x(Li) = 0.98) P3 The /3 phase is consistently the low- 
est energy phase because Z, the most electronegative el- 
ement, is nn to both Mn and Li. 

The bonding properties can be analyzed using the 
states at the r-point of the BZ because these states most 
clearly reflect their atomic origins. Fig. [2] shows the den- 
sity of states (DOS) for the a and j3 phases of LiMnP 
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TABLE II. The optimized lattice constants a opt , total ener- 
gies i?Tot, and magnetic moments M for three lithiated half- 
Heusler alloys. 



TABLE III. Summary of the lattice constants a,h m , total en- 
ergies Exot, magnetic moments per unit cell M, and energy 
gaps A near the Fermi energy Ef (with spin channel) of the 
three phases where the magnetic moment is largest, or integer. 



with markers indicating the energy and atomic contri- 
bution to the states at the T-point. Here, we choose to 
only discuss the bonding properties of LiMnP since, at 
the optimized lattice constant, not many properties are 
significantly altered for the other two compounds. Also, 
the j3 and 7 DOSs are very similar — both have Mn and P 
as nn and the position of the Li atom does not drastically 
alter the DOS — so the 7 phase is not included. The f3 
and 7 phases form a large bonding/antibonding gap, re- 
sulting from the hybridization of the so-called t2 9 states 
(the triply degenerate d xy , d yz , d zx states) of Mn and 
p states (pz, p y , p z ) of P, because Mn and P are nn in 
the tetrahedral environment. The a phase is significantly 
different than the /3 and 7 phases. This structure does 
not form a large bonding/antibonding gap because the 
Mn and P are second neighbors in a cubic environment. 



B. Can the three half-Heusler alloys be 
ferromagnetic or half-metallic? 

We are interested in adjusting the lattice constant to 
force the values of the magnetic moments to be as large 
as possible. In practice, this is possible because many de- 
vices are thin films and can be grown to match the lattice 
constants on selected substrates. We can use our results 
to predict the substrates that match the lattice constants 
of our compounds with large magnetic moments, but this 
will not always be possible for large lattice constants or 
energies far from equilibrium. Tetragonalization can also 
occur during growth, but it wil l no t necessarily destroy 
the half-metallicityP31 In Table III we present detailed 



ferromagnetic properties of the three compounds in the 
primitive cells of the three types of structures. 

If we only consider the primitive unit cell, then ev- 
ery alloy in Table |III| with an integer magnetic moment 
and a gap is a half-metal. These include: a-LiMnN, a- 
LiMnP, /3-LiMnSi /3-LiMnP, and 7-LiMnSi. The lattice 
constant can be increased slightly without destroying the 
half-metalicity since there is a range of lattice constants 
where Ep falls within the gap of the semiconducting gap. 



The magnetic moments of the half-metals agree with the 
ionic model plus the contribution of the Li electron to the 
moment of Mn. The remaining elements are ferromag- 
nets and any change in the lattice constant will decrease 
the moment. Some of the important features of the indi- 
vidual half-Heusler alloys are discussed below. 

Fig. [3] shows the band structures of 7-LiMnSi at two 
lattice constants: the optimized lattice constant and the 
half- metallic lattice constant. The lowest energy band is 
formed by the s-state of the Z atom and is not shown. 
The next bands, shown in the figures, are the t^g states 
that split from triply degenerate states at the T-point to 
doubly (upper) and singly (lower) degenerate states as 
k moves towards X. The next higher energy states are 
originally from e g states of Mn. Since their lobes point 
away from nn Si and toward the sn of the Mn.P they do 
not strongly interact with any other states so they are 
called the non-bonding states. The half-metallic gap, in 
Fig. j3^b) , is formed between the e g states and the doubly 
degenerate states that split off of the t2 g states. Using the 
band structure, we can understand why this half-Heusler 
alloy is not a half metal at the optimized lattice constant. 
The e g states are insensitive to the separation between 
nn, but the smaller optimized lattice constant causes Ep 
to shift up and intersect with the e g states resulting in 
the disappearance of the half-metalicity. 

The reason that a-LiMnSi is not a half-metal can be 
attributed to the fact that Si forms nn pair with Li and 
sn pair with Mn. The qualitative details are: (i) Ideally, 
the Li electron can form a covalent bond with Si because 
the electronegativity of Si is not as large as either P or 
N. (ii) The sn configuration between Si and Mn favors 
the p-states of Si to form bonds with three Mn electrons. 
The strengths of all these bonds are not large enough 
to form a gap. On the other hand, at the TM site, the 
remaining four d-electrons under Hund's rule align their 
spins and give 4^b- 

The pnictides in the 7 structure do not have integer 
magnetic moments and show no gap in either spin chan- 
nel near Ep. According to our criteria,^ they are not 
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FIG. 2. Spin polarized DOS for LiMnP in (a) a and (b) j3 phases. Labels indicate the peaks and the fraction of atomic states 
per band (insets), at the T-point, occupying p-states (p x , p y and p z ) around P or Li, and e 9 (d z 2 and d x 2_ y 2) and t2 g {d xy , 
d yz , and d zx ) states of Mn. For clarity, the origins of other peaks in the DOS are not included. 




FIG. 3. (Color) F-X Band structure for 7-LiMnSi at (a) the optimized lattice constant a — 5.715 A, and (b) the half-metallic 
lattice constant a = 6.749 A. The triply degenerate states, t2 9 , and the doubly degenerate states, e 9 , near Ep are labeled. 



6 



half metals. To see the effect of Li in pnictide LiMnZ, 
we present in Fig. [4] the DOS for (a) zinc blende MnP 
and (b) 7-LiMnZ along with the size of the bonding- 
antibonding gap (the energy difference between the top of 
the valence band and the bottom of the conduction band, 
in the minority spin channel, at the T-point). The Li 
electron is not completely transferred to the Mn because 
of the presence of the highly electronegative pnictogen. 
Instead, the electron from Li causes a reduction in the 
strength of the d-p hybridization between the pnictogen 
and the Mn under the tetrahedral environment. With the 
reduction in the hybridization, the bonding-antibonding 
gap shrinks and the bottom of the conduction band in the 
minority spin channel of 7-LiMnP, or 7-LiMnN, becomes 
occupied. 

The pnictides exhibit half metallic properties in the 
a structure at lattice constants larger than their respec- 
tive optimized values. The gap for a-LiMnP is 0.376 eV 
while it is 0.696 eV for a-LiMnN. These values reflect the 
strength of the electronegativity of P with respect to N 
and the differing lattice constants. The gaps are formed 
by the t2 S -states of Mn and the p-states of the pnicto- 
gen, but are much smaller than the gaps in the /3 and 7 
phases. The DOS for a-LiMnP is provided in Fig. [5] 

Based on the fact that the density of states of 7-LiMnP, 
given in Fig. |4]^b), shows E F slightly above the conduc- 
tion band edge in the minority spin channel, we decided 
to find new half metals by removing one of the Li from 
the conventional cell. The symmetry of the half-Heusler 
phases allows the removal of any Li atom in the unit 
cell. The resultant sample is labeled as 7~Li3Mn 4 P4 in 
the conventional cell or as 7-Lio.75MnP as an alloy. After 
the ionic relaxation resulting from the removal of a Li, we 
found a half metal at the lattice constant of 6.57 A. The 
spin-polarized density of states, shown in Fig. [6j clearly 
indicates the semiconducting feature of a half metal — Ep 
is in the minority spin channel gap. The corresponding 
magnetic moment is 19.000 /ie /unit cell. This result sug- 
gests that Li is giving an electron to Mn in the 7 struc- 
ture: the missing Li reduces the magnetic moment from 
20 to 19 /iB/unit cell. A detailed table of the remaining 
phases, with one Li missing, at the half metallic or max- 
imum magnetic moment lattice constant is presented in 
Table IV In the a phase of Li3Mn 4 N4, no maximum or 



Phase Compound ai lm (A) Erot (eV) M (/xb/u.c.) A 



half metallic moment occurred for any lattice parameter 
between 4.00 A and 8.00 A so we did not include it in Ta- 
ble |IV| Outside this range, the total energy is too high 
to be physical. 

Many of the other materials in Table |IV| obey the 
simple behavior of reducing the magnetic moment of 
Li 4 Mn4Z4 by one, but there are a few exceptions. We 
shall dicuss the exception in the a phase in the next sec- 
tion, but 7-Li 3 Mn 4 N4 in the conventional cell has a mag- 
netic moment of only 17 /j,b /unit cell. Also, the lattice 
constant of 7-LiaMn4N4 is small corresponding to the 
large electronegativity of N. The confined spin results in 
the alignment of the spins, similar to what we previously 
found with MnCp^ and a lower magnetic moment. 
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Li 3 Mn 4 Si 4 
Li 3 Mn 4 P 4 


5.83 
5.91 


-59.695 
-64.133 


14.078 
13.000 


Metal 
0.094 | 






5.64 


-70.593 


19.000 


1 905 




Li 3 Mn 4 Si 4 


6.25 


-59.434 


15.000 


0.903 | 




Li 3 Mn 4 P 4 


6.24 


-64.418 


19.000 


1.747 X 


7 


Li 3 Mn 4 N 4 
Li 3 Mn 4 Si 4 
Li 3 Mn 4 P 4 


5.33 
6.40 
6.57 


-70.746 
-57.585 
-60.169 


17.000 
15.009 
19.000 


1.127 i. 
Metal 
0.982 | 



TABLE IV. Summary of the lattice constant Ohm, total en- 
ergy i/rot, magnetic moment per formula-unit M, and en- 
ergy gap A (or "Metal" if there is no gap) of Li 3 Mn 4 Z 4 in 
the three phases. a-Li 3 Mn 4 N 4 is not included because it re- 
mained metallic between a = 4.00 A and 8.00 A, and we did 
not observe half-metallicity or a maximum moment. 



C. Should the Slater-Pauling curve, originally 
formulated for the TM-based half-Heusler alloys, be 
altered for the lithiated alloys? 

The lithiated half Heusler alloys either have 13 or 12 
valence electrons per unit cell. Since N t is smaller in these 
alloys than in the TM-related compounds, the Slater- 
Pauling curve should definitely be modified. The band 
structure of 7-LiMnSi in the primitive cell, shown in Fig. 
|4^b), shows four bands in the semiconducting channel, 
N± = 4. Eq. (TH) must be modified: 



M = (Nt — 8) /xs/f. u. 



(2) 



For LiMnSi, N t = 12 and V; = 4, so M is 4> B /f.u. 
Similarly, M = 5 /zb /primitive cell for the pnictides. The 
calculated moments agree well to the predicted results of 
Eq. (2j. 

Most of the LiaMn 4 Z 4 compounds agree with the 
Slater-Pauling curve if we subtract one to account for 
the missing electron. The conventional cells have ei- 
ther Nt = 4 x 13 — 1 = 51, for the pnictides, or 
N — 4x 12 — 1 = 47, for Li3Mn 4 Si 4 , valence electrons per 
unit cell. According to our modified rule, the moment for 
the conventional cell should be 



M= (JV t -32)/i B /f.u. 



(3) 



which gives 15/iB/f-u. for Li3Mn 4 Si 4 and 19/XB/f. u. for 
the pnictides. This expression agrees with our results 
for the P and 7 phases, however, as seen in Table |IV| 
Eq. ^ does not agree with results for the a structure. 
The disagreement is a consequence of the relaxation of 
the 4 Mn and 4 P atoms surrounding the vacancy left 
by removing the Li. The Mn and Z atoms move toward 
the vacancy, a distance of 0.96 A and 0.24 A, respectively, 
form bonds with each other and the Mn d-electrons pair 
up, reducing the overall moment. This is the same spin 
confinement effect found in MnO^and 7-Li 3 Mn 4 N 4 . As 
a result of the different configuration of up and down spin 
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FIG. 4. Density of states for (a) zinc blende MnP and (b) 7-LiMnP conventional cells. Points (i)-(iv) indicate the energy and 
atomic contribution, per band, to the bonding and antibonding states (inset), at the P-point, in the minority spin channel. 
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FIG. 5. The DOS of a-LiMnP showing the bonding-antibonding gap between points (iv) and (v) and the contribution of the 
projected states — p, t 2g and e g — around each atom (inset). 
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FIG. 6. DOS for LisMn4P4. The inset specifies the atomic character of the peaks near the points that make up the top of the 
valence band (i) and the bottom of the conduction band (ii) at the T-point in the minority spin channel. 



Phase 


Compound 


a (A) Bafm (eV) 


E FM (eV) 


AE (eV) 




Li 2 Mn 2 N 2 


5.491 


-33.115 


-33.078 


-0.037 


a 


Li 2 Mn 2 Si 2 


5.975 


-30.208 


-30.308 


+0.100 




Li 2 Mn 2 P 2 


6.525 


-29.963 


-29.937 


-0.026 




Li 2 Mn 2 N 2 


5.641 


-37.202 


-36.910 


-0.292 


P 


Li 2 Mn 2 Si 2 


6.590 


-30.046 


-30.393 


+0.347 




Li 2 Mn 2 P 2 


6.809 


-31.858 


-31.735 


-0.123 




Li 2 Mn 2 N 2 


5.571 


-34.481 


-34.179 


-0.302 


7 


Li 2 Mn 2 Si 2 


6.749 


-28.594 


-28.918 


+0.324 




Li 2 Mn 2 P 2 


7.071 


-29.130 


-29.043 


-0.087 



TABLE V. Summary of lattice constant a, total energies of 
antiferromagnetic -Eafm and ferromagnetic -Efm orderings, 
and the energy difference AE between AFM and FM order- 
ings of the three compounds. Negative energy differences in- 
dicate that the AFM phase is energetically preferred. 



electrons, the band structure is different and the premise 
behind Eq. ([3| is flawed. In the j5 and 7 structures, since 
the Mn and Z atoms already form strong bonds, remov- 
ing the Li causes the other atoms to move only 0.04 A, 
essentially maintaining the same nn bond distance and 
not confining the spins. 



D. Can the three alloys be fully compensated (FC) 
half-metals or antiferromagnets? 

There is only one magnetic transition metal element 
in the formula unit, so it is unlikely to find the FC state 
in these materials. Equations ^ and ^ also suggest 
we will not find the FC state since there are more va- 
lence electrons (N t = 12 or 13) than are likely to align 
antiparallel (2ATj, = 8). Instead, we focus on the antifer- 
romagnetic cases. In order to investigate this phase, we 
use the aforementioned tetragonal cell. In Table |VJ we 
present the results for the three half Heusler alloys in the 
antiferromagnetic phase and their energy with respect to 
the ferromagnetic case. 

From our calculations, all three phases of LiMnSi can 
be ferromagnets while LiMnN and LiMnP favor he an- 
tiferromatic phase at the half-metallic lattice constant, 
much like LiMnAsP 



E. How does Li play its role in the two structures? 

We have argued that Li easily gives up its valence elec- 
tron to its nearest neighbor. To substantiate this argu- 
ment, we compare 7-LiMnP to MnP in the ZB structure. 
The latter is predicted to be a half metal in which the 
d-p hybridization and the non-bonding states are easily 
identified.^ I n Fig. Qa) and (b) , we compared the states 
around Ep at the T-point of the two spin channels for 
LiMnP and MnP at the same lattice constant of 6.57 A. 
The argument that the Li is giving up its electron is re- 
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FIG. 7. (Color) Charge density difference of 7-LiMnP-ZB 
MnP at the same lattice constant. Colored circles indicate 
locations of the atoms. The (i) and (ii) labels indicate where 
the Li electron bonds with Mn and P, respectively. Label (iii) 
shows the effect of the d-p hybridized bond weakening due to 
the presence of the Li electron. 



fleeted in the character of the electron around Li. The 
Li s-state is partially promoted to the p-state, at points 
(iii) and (iv) of Fig. Efb), and contributes to the d-p 
mixed states. Furthermore, this fact is demonstrated by 
examining the charge density difference, LiMnP— MnP, 
shown in Fig. [7j This plot shows that the Li electron 
redistributes its charge towards Mn (i) and P (ii), away 
from Li, and causes the d-p hybridized bond between Mn 
and P to weaken, thereby shifting the bond charge away 
from P and towards Mn (iii). 



F. Can the lithiated alloys in the /3-phase be more 
stable than the corresponding compounds in the 
zinc blende structure? 

In the ZB structure, a weak restoring force constant in 
response to a shear stress is demonstrated by the soften- 
ing of the transverse acoustic (TA) branch of the phonon 
energy at the zone boundary along the [100] direction, 
the X pointP^We carried out response-function phonon 
calculations for /3-LiMnP and ZB MnP: the results are 
shown in Fig. [§} /3-LiMnP has increased stability over 
ZB MnP along the [100] direction and shows stability in 
the [110] direction. MnP is unstable in the [110] direc- 
tion. 
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FIG. 8. The phonon bands on K-T-X of (a) /3-LiMnP and (b) zinc blende MnP. 



IV. SUMMARY 

We investigated three Li-based half Heusler alloys, in- 
volving one TM element per formula, for the possibil- 
ity of being ferromagnetic with large magnetic moments. 
Three different arrangements of the atoms, denoted a, f3 
and 7, appeared in the literature. The 7 structure can 
be viewed as an intercalation of Li between the other two 
elements. We expect that many properties differ from 
the TM-related half Heulser alloys involving two TM ele- 
ments. We compare the bonding properties of LiMnP, a 
prototype, to MnP, a half metal in the zinc blende struc- 
ture. Using the tetragonal cell, a-LiMnP, ev-LiMnN and 
/3-LiMnP are antiferromagnetic while LiMnSi is a ferro- 
magnetic half-metal in the j3 and 7 phases, all at larger 
lattice constants than their respective optimized lattice 
constant. The reasons for the above facts are provided. 
The magnetic moments for the ferromagnetic alloys are 
larger than 3 /is /formula unit. The conventional cell has 
been used to find Li 3 Mn4P4 and Li 3 Mn4N 4 are half met- 
als along with Li3Mn4Si4 in the /3 phase. A new curve 



predicting the magnetic moments in this class of half 
Heusler alloys is proposed and predicts well for the three 
alloys in the ferromagnetic phase. The stability against 
shear stress as compared to the simple ZB structure is 
demonstrated by calculating the phonon spectrum, along 
the [100] direction of the Brillouin zone, using the phonon 
response function scheme. The added Li in the half- 
Heusler alloy stabilized the structure with respect to the 
zinc blende structure. In addition, we show that MnP 
is unstable along the [110] direction. Hopefully, these re- 
sults will facilitate the search of new half metals involving 
alkali metals. 
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